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Abstract 

A  model  dealing  with  the  anode  catalyst  contamination  induced  by  fuel  impurities  has  been  developed.  This  model  can  be  used  to  describe 
the  transient  and  steady-state  performance  losses.  Several  characteristics  such  as  performance  loss,  contamination  transient  time  constant  and 
recovery  process  have  also  been  introduced  into  the  model.  The  obtained  equations  can  be  used  to  simulate  and  estimate  the  chemical  and 
electrochemical  reaction  rate  constants,  and  make  some  prediction  about  the  severity  of  the  contamination  and  the  performance  recoverability. 
Crown  Copyright  ©  2005  Published  by  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  contamination  caused  by  the  impurities  in  the 
feed  stream  is  one  of  the  important  issues  in  fuel  cell  devel¬ 
opment  and  operation  [1-24].  The  major  impurities  in  the 
hydrogen  stream  are  CO,  F^S,  NH3,  organic  sulfur-carbon 
and  carbon-hydrogen  compounds,  etc.,  produced  in  the  hy¬ 
drogen  production  process  (for  example,  natural  gas  re¬ 
forming).  These  impurities  (or  contaminants)  degrade  the 
cell  performance  and  sometimes  cause  permanent  membrane 
electrode  assembly  (MEA)  damage.  The  two  effects  can  be 
used  to  describe  the  contamination  process:  (1)  kinetic  ef¬ 
fect  (poisoning  of  the  electrode  catalyst)  and  (2)  conductivity 
effect  (increase  in  the  solid  electrolyte  resistance  including 
those  of  membrane  and  catalyst  layer  ionomer).  It  is  well 
known  and  widely  documented  that  the  major  impact  of 
CO  contaminant  on  the  MEA  performance  is  the  kinetic  ef¬ 
fect  that  is  the  poisoning  of  the  anode  catalyst  [1,25-27]. 
For  H2S  contamination,  the  impact  is  also  kinetic  related 
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and  has  also  been  recognized  as  a  severe  depression  in 
the  fuel  cell  performance  for  many  years  [1,3,12,14,18-24]. 
NH3,  as  one  of  the  fuel  impurities,  has  also  been  investi¬ 
gated  in  fuel  cell  operation  [2,3]  and  the  conductivity  ef¬ 
fect  has  been  found  to  be  the  major  cause  for  the  fuel  cell 
degradation. 

It  has  been  found  that  the  performance  loss  is  related  to 
the  contaminant  concentration  in  the  feed  stream  [2,3],  and 
the  fuel  cell  operating  current  density  (or  cell  voltage).  In 
general,  the  higher  the  contaminant  level  or  current  density, 
the  faster  and  deeper  the  depression  is.  For  the  purpose  of 
fundamental  understanding  and  application  development,  a 
model  is  necessary  to  describe  theoretically  the  contamina¬ 
tion  process.  For  the  CO  poisoning  process,  several  model 
papers  have  been  published,  which  mainly  focused  on  the 
steady-state  poisoning  effect  [26,28-30].  For  a  transient  pro¬ 
cess,  Bhatia  and  Wang  has  reported  a  model  and  validated  it 
by  experimental  results  [30] .  A  general  and  complete  model 
for  the  contamination  process  as  a  function  of  contaminant 
level  and  current  density,  including  transient,  steady-state  and 
recovery  processes,  is  expected  to  be  very  useful  in  fuel  cell 
research  and  development. 
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Nomenclature 

a  as  defined  in  the  text  for  the  recovery  process 

(mol3  cm-6  h-3) 

b  as  defined  in  the  text  for  the  recovery  process 

(mol3  cm-6  h-3) 

A  as  defined  in  the  text 

Ao  as  defined  in  the  text 

B  as  defined  in  the  text 

B\  sum  of  the  rate  constants  (mol  cm-2  h-1) 

#2  backward  rate  for  reaction  (7)  (mol  cm-2  h- 1 ) 
B\  sum  of  the  rate  constants  at  steady- state  anode 

overpotential  (mol  cm-2  h-1) 

B I  backward  rate  for  reaction  (7)  at  steady-state 

anode  overpotential  (mol  cm-2  h-1) 

C  as  defined  in  the  text 

Ch2  concentration  of  hydrogen  in  the  electrolyte 

(Nafion  ionomer  catalyst  layer)  (mol  cm-3) 
Co  concentration  of  air  oxygen  in  the  electrolyte 
(Nafion  ionomer  catalyst  layer)  (mol  cm-3) 
CH+  proton  concentration  in  the  catalyst  layer 
(mol  cm-3) 

Cp  contaminant  concentration  (ppm) 

Cp'  contaminant  reaction  product  concentration 

(ppm) 

£a  anode  potential  at  7a  ^  0  and  Cp  ^  0  (V) 

E^  anode  potential  at  7a  =  0  and  Cp  =  0  (V) 

£0a  standard  anode  potential  (V) 

Ec  cathode  potential  at  7a  ^  0  and  Cp  ^  0  (V) 

E [?  cathode  potential  at  7a  =  0  and  Cp  =  0  (V) 

F  Faraday’s  constant  (96,487  A  s  mol-1) 

F\  sum  of  the  forward  rate  constants 

(mol  cm-2  h-1) 

F2  forward  rate  for  reaction  (7)  (mol  cm-2  h- 1 ) 
Ff  sum  of  the  forward  rate  constants  at  steady- 
state  anode  overpotential  (mol  cm-2  h-1) 
forward  rate  for  reaction  (7)  at  steady- state  an¬ 
ode  overpotential  (mol  cm-2  h-1) 

F[~t0  sum  of  the  forward  rate  constants  after  the  con¬ 

taminant  is  cut-off  (mol  cm-2  h-1) 
exchange  current  density  for  hydrogen  oxida¬ 
tion  on  the  platinum  surface  in  the  absence  of 
contaminant  (A  cm-2) 

/q  exchange  current  density  for  oxygen  reduction 
on  the  platinum  surface  in  the  absence  of  con¬ 
taminant  (A  cm-2) 

4  fuel  cell  anode  current  density  (A  cm-2) 

4  fuel  cell  cathode  current  density  (A  cm-2) 

/Ceii  fuel  cell  operating  current  density  (A  cm-2) 
k\f  reaction  (1)  forward  reaction  rate  constant 
(cmh-1) 

km  reaction  (1)  backward  reaction  rate  constant 
(mol  cm-2  h-1) 


4>f  reaction  (2)  forward  reaction  rate  constant 
(mol  cm-2  h-1) 

4>b  reaction  (2)  backward  reaction  rate  constant 
(mol  cm-2  h-1) 

&3f  reaction  (3)  forward  reaction  rate  constant 
(mol  cm-2  h-1) 

&3b  reaction  (3)  backward  reaction  rate  constant 
(cmh-1) 

&4f  reaction  (4)  forward  reaction  rate  constant 
(mol  cm-2  ppm- 1  h- 1 ) 

&4b  reaction  (4)  backward  reaction  rate  constant 
(mol  cm-2  h-1) 

&5f  reaction  (5)  forward  reaction  rate  constant 
(mol  cm-2  ppm- 1  h- 1 ) 

4)b  reaction  (5)  backward  reaction  rate  constant 
(cmh-1) 

&6f  reaction  (6)  forward  reaction  rate  constant 
(mol  cm-2  ppm- 1  h- 1 ) 

kbb  reaction  (6)  backward  reaction  rate  constant 
(cmh-1) 

&7f  reaction  (7)  forward  reaction  rate  constant 
(mol  cm-2  h-1) 

^7b  reaction  (7)  backward  reaction  rate  constant 
(mol1-q  cm3q-2  h-1) 

&8f  reaction  (8)  forward  reaction  rate  constant 
(mol  cm-2  h-1) 

&8b  reaction  (8)  backward  reaction  rate  constant 
(mol-q  cm3q+ 1  h-1) 

m  water  stoichiometry  in  reactions  (7)  and  (8) 
nao  electron  transfer  coefficient  for  overall  cathode 
oxygen  reduction 

nn  electron  transfer  number  for  overall  hydrogen 
oxidation 

no  electron  transfer  number  for  overall  cathode 
oxygen  reduction 

7*3,5, 7-8  electron  transfer  number  for  individual  elec¬ 
trochemical  half-reaction 
P  symbol  for  contaminant 

P;  symbol  for  the  product  of  P  electrochemical 

oxidation 

PL%  fuel  cell  performance  loss  in  the  presence  of 
contaminant 

R  Gas  constant  (8.314  JK-1  mole-1) 

q  proton  and  electron  numbers  in  reactions  (7) 

and  (8) 

Ro  fuel  cell  internal  resistance  (electrolyte  (mem¬ 

brane)  resistance)  (Q  cm2) 

to  lifetime  moment  at  which  the  contaminant 

source  is  cut-off  (h) 

T  temperature  (K) 

Vceii  fuel  cell  voltage  (V) 
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C  =0 

Kell  steady  fuel  cell  voltage  in  the  absence  of  con¬ 
taminant  (V) 

Vcsell  fuel  cell  steady- state  cell  voltage  in  the  pres¬ 
ence  of  contaminant  (V) 

V°  fuel  cell  open  circuit  voltage  (V) 

Z  as  defined  in  the  text  (mol  cm-2  h_  1 ) 

Z1-1  o  as  defined  in  the  text  (mol  cm  2  h  1 ) 

Greek  symbols 

<^3,5, 7-8  electron  transfer  coefficient  for  individual 
electrochemical  half-reaction 
ao  electron  transfer  coefficient  for  fuel  cell  cath¬ 
ode  oxygen  reduction 

ya  real  active  surface-to-electrode  geometric 
surface  ratio  of  the  anode  catalyst  layer 
(cm2  cm-2) 

y c  real  active  surface-to-electrode  geometric  sur¬ 

face  ratio  of  the  cathode  catalyst  layer 
(cm2  cm-2) 

r pt  Pt  surface  site  concentration  (mol  cm-2) 
r pt_H  Pt  surface  site  concentration  (mol  cm-2) 
7pt-H2  Pt  surface  site  concentration  (mol cm-2) 
r pt_p  Pt  surface  site  concentration  (mol  cm-2) 
7"pt_p'  Pt  surface  site  concentration  (mol  cm-2) 

/pt  total  surface  site  concentration  (mol  cm-2) 

?7a  anode  overpotential  (V) 

77c  cathode  overpotential  (V) 

0pt  surface  coverage  of  platinum  unoccupied  site 

0pt  surface  coverage  of  platinum  unoccupied  site 

at  equilibrium  anode  potential 
$Pt_H  surface  coverage  of  hydrogen  atom  on  the  plat¬ 

inum 

0pt_H  hydrogen  atom  surface  coverage  on  platinum 
at  equilibrium  anode  potential 
0Pt_ h2  hydrogen  molecule  surface  coverage  on  the 
platinum 

0pt_p  contaminant  surface  coverage  on  the  platinum 
0pt_P'  contaminant  reaction  product  surface  coverage 

on  the  platinum 

0pf_!p  contaminant  surface  coverage  on  the  platinum 
after  the  contaminant  source  is  cut-off 
0ptL.p/  contaminant  reaction  product  surface  coverage 

on  the  platinum  after  the  contaminant  source 
is  cut-off 

Cp  contaminant  surface  coverage  on  the  platinum 
at  the  moment  of  t  =  to 

6*-%,  contaminant  reaction  product  surface  coverage 

on  the  platinum  at  the  moment  of  t  =  to 
r  contamination  transient  time  constant  (h) 


In  this  paper,  only  the  kinetic  effects  will  be  considered  in 
the  modeling  process.  Thus,  a  model  dealing  with  the  gen¬ 
eral  case  of  anode  catalyst  contamination  by  fuel  impurities 
has  been  developed.  This  model  can  be  used  to  describe  the 
transient  and  steady-state  performance  loss  processes.  Sev¬ 
eral  characteristic  parameters  such  as  performance  loss,  con¬ 
tamination  transient  time  constant,  and  recovery  period  have 
also  been  introduced  to  describe  the  contamination/recovery 
processes. 

2.  Model  description 

2.7.  Proposed  anode  chemical  reactions  in  the  presence 
of  contaminant  P 

In  order  to  simulate  the  degradation  behavior  of  the  cell 
performance  in  the  presence  of  a  stream  contaminant,  and 
determine  quantitatively  the  relationship  of  contamination 
surface  coverage  as  a  function  of  the  contaminant  concentra¬ 
tion,  cell  current  density  and  lifetime,  several  possible  chem¬ 
ical  and  electrochemical  reactions  with  their  correspond¬ 
ing  reaction  constants  have  been  assumed,  as  in  reactions 
(I)— (VIII): 

f  f 

Pt  +  H2  <=*  Pt-H2  (I) 

klb 

kif 

Pt-H2  +  Pt  ^  2Pt-H  (II) 

klb 

%  &xp(a3n3Fqa/RT) 

Pt-H  +±  Pt  +  H+  +  e-  (III) 

&3b  exp(— (1— ar3)ft3F77a/ RT) 

&4f 

P  +  Pt  <=►  Pt-P  (IV) 

&4b 

k5f  exp(a5n5Fqa/RT) 

P  +  Pt-H  =:  Pt-P  +  H++e"  (V) 

k5h  exp(-(l-a5)n5Fqa/RT) 

P  +  Pt-H2  <=>  Pt-P  +  H2  (VI) 

hb 


Pt-P  +  mH2Q 


&7f  exp(or/qFr)a/RT ) 


Pt-Pr  +  gH+  +  qt 


k7h  exp(-(l-or/)qFria/ RT) 


(VII) 


%  expia^q Frill / RT) 

Pt-P  +  mH2Q  Pt  +  P  +  gH+  +  qt 


h b  exp(— (1  -as)qFqa/RT) 


(VIII) 


Reactions  (I)— (III)  are  the  mechanism  for  hydrogen  oxi¬ 
dation,  which  has  been  investigated  and  reported  for  many 
years  [25,26,31-44].  Reaction  (I)  is  the  adsorption  of  H2  on 
the  platinum  surface,  which  is  a  fast  reaction.  Reaction  (II) 
is  the  slow  dissociative  chemical  adsorption  of  adsorbed  H2, 
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which  has  been  considered  as  the  rate-determining  step  for 
H2  oxidation  [25,26]  and  reaction  (III)  is  the  fast  electro¬ 
chemical  oxidation  of  the  dissociated  hydrogen  atom.  Re¬ 
actions  (IV)-(VI)  are  those  surface  adsorption  and  surface 
electrochemical  reaction  of  the  poisoning  species  (marked  as 
P).  Note  that  reactions  (V)  and  (VI)  proposed  here  are  used 
to  describe  the  interaction  between  the  contaminant  species 
P  and  the  atomic  and  molecular  hydrogen  occupied  Pt  sur¬ 
face.  It  is  worthwhile  to  point  out  that  in  some  situations, 
one  contaminant  P  molecule  could  react  with  several  surface 
Pt  atoms  to  form  surface  species  such  as  Ptn-P  (i.e.,  n  >  1). 
However,  in  order  to  make  the  modeling  process  simpler, 
only  Pt-P  is  considered  as  the  surface  poisoning  species  in 
this  paper.  Reaction  (VII)  or  (VIII)  is  the  oxidation  of  ad¬ 
sorbed  contaminant  on  the  Pt  surface,  resulting  in  a  prod¬ 
uct  that  could  be  either  a  surface  adsorbed  species  (Pt-P') 
or  a  soluble  species  (Pr).  In  the  case  that  P  is  not  electro- 
chemically  reactive,  the  corresponding  reactions  (VII)  and 
(VIII)  would  be  removed  from  the  proposed  mechanism. 
Each  reaction  has  its  own  rate  constants  (%  for  forward  and 
k/b  for  backward).  For  those  electrochemical  reactions  (re¬ 
actions  (III),  (V),  (VII)  and  (VIII)),  the  rate  constants  are 
electrode  overpotential  dependent,  and  written  according  to 
the  Butler- Volmer  equation.  Where  c^,  rii  and  q  are  electron 
transfer  coefficient,  electron  transfer  number  and  the  pro¬ 
ton/electron  numbers  for  the  corresponding  electrochemical 
reaction. 

2.2.  Surface  coverage  and  their  expressions  as  a 
function  of  anode  overpotential 

The  anode  overpotential  is  the  difference  between  two 
stages  of  anode  potential,  i.e.,  the  anode  equilibrium  state 
potential  (F^)  at  which  the  anode  net  current  density  (If) 
and  the  contaminant  concentration  (Cp)  are  both  equal  to 
zero,  and  the  anode  potential  (Fa)  at  which  both  /a  and  Cp 
are  not  equal  to  zero,  that  is,  r)a  =  Ea  —  E®.  In  order  to  ob¬ 
tain  the  time  dependence  of  contaminant  surface  coverage, 
the  interrelationship  among  various  kinds  of  surface  cov¬ 
erage  and  the  anode  overpotential  have  to  be  obtained.  As 
in  the  reaction  mechanism  proposed  above,  there  are  five 
kinds  of  surface  sites,  which  are  Pt,  Pt-H2,  Pt-H,  Pt-P  and 
Pt-P2  Their  corresponding  surface  concentration  can  be  ex¬ 
pressed  as  r pt,  Cpt_H,  Tpt— h2  ,  Cpt_p,  and  FPt_P'  with  a  unit  of 
mol  cm-2 .  The  sum  of  these  surface  concentrations  can  be  de¬ 
fined  as  the  total  available  Pt  surface  concentration,  Fpj  with 
a  unit  of  mol  cm-2 .  The  surface  coverage  for  each  surface 
site  can  be  expressed  as  0Pt  (=FPt/FpTt),  0Pt_H(=  F^-h/F^), 

$Pt— h2(=  fpt-H2/^Pt)’  ^Pt-p(=  fpt-p/fpi)  and  $Pt-P'(= 

FPt_p//FpJ),  respectively.  Fig.  1  shows  the  schematic  ex¬ 
pression  of  a  platinum  surface  with  different  surface 
sites. 

The  equation  (1)  can  be  used  to  describe  the  inter¬ 
relationship  among  the  surface  coverage  terms: 


Fig.  1 .  Schematic  diagram  of  a  platinum  surface  with  adsorption  of  molecu¬ 
lar  hydrogen,  atomic  hydrogen,  contaminant  and  contaminant  reaction  prod¬ 
uct. 


At  anode  equilibrium  potential  in  the  absence  of  contam¬ 
inant,  equation  (2)  can  be  used  to  describe  the  surface 
coverage: 

0pt  +  #pt_H  +  #Pt-H2  =  1 

Since  reaction  (II)  is  the  rate-determining  step,  it  can  be 
assumed  that  reaction  (I)  is  always  at  its  equilibrium  state 
whether  there  is  a  contaminant  present  or  not.  From  reaction 
(I),  the  relationship  between  0 pt  and  0pt-H2  can  be  establi¬ 
shed: 


where  Ch2  is  the  average  concentration  of  hydrogen  in  the 
vicinity  of  the  catalyst  layer,  expressed  either  by  gas  concen¬ 
tration  or  the  wet  concentration  in  the  electrolyte  (ionomer 
matrix  layer  in  the  catalyst  layer,  mol  cm-3). 

If  reaction  (III)  is  very  fast,  the  anode  potential  would 
follow  the  Nernst  behavior  even  in  the  presence  of  poisoning 
species  [25,26]: 


Fa  =  F0a  + 


RT 

~E 


In 


RT 

H - In 

F 


(On  CH+\ 

V  0Pt— H  ) 


where  F°a  is  the  standard  anode  potential  of  the  hydrogen 
redox  reaction  and  CH+  is  the  concentration  of  proton.  At 
zero  current  density  (reaction  (III)  at  its  equilibrium  state), 
and  if  the  contaminant  species  P  is  absent,  the  equilibrium 
electrode  potential  would  be: 


A  A  RT 

E0  =  £0a  - ln 

a  F 


RT 

H - In 

F 


$Pt  +  ^Pt— H  +  $Pt— H2  +  0Pt-P  +  $Pt— P'  —  1 


(5) 
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where  0pt  and  0pt_H  are  the  surface  coverage  at  equilibrium 
anode  potential  in  the  absence  of  the  contaminant.  Com¬ 
bining  equations  (l)-(5),  the  surface  coverage  of  hydrogen 
molecules,  hydrogen  atoms  and  unoccupied  Pt  sites  can  be 
expressed  as  a  function  of  anode  overpotential  r)a  and  con¬ 
taminant  coverage  0pt_p  and  0pt_p' : 


(9Pl  =  A(\  -  0Pt_p  -  0Pt_p') 


(6) 


where 


A  = 


t^pt 


O' 


o 

Pt— H 


exp 


Jjh 

RT 


1  + 


^Pt 


o 

Pt— H 


O' 


1  +  ^Ch2  )  exp 


Fjh 

RT 


(7) 


(8) 


$Pt-H2  =  B(  1  —  #Pt-p  —  Opi-Pf) 
where  B  =  |^Ch2A,  and 
0pt-H  =  C(1  —  0pt-p  -  #pt-p') 

where  C  =  1  —  ^1  +  A. 

The  surface  coverage  ratio  in  equation  (6),  0pt/0pt_H,  is 
determined  by  the  equilibrium  anode  potential  and  proton 
concentration  as  shown  in  equation  (5). 


2.3.  Contaminant  surface  coverage  as  a  function  of  r\a 
and  lifetime  (t) 

From  reaction  (IV)  to  (VII),  the  increase  in  contaminant 
surface  coverage  with  time  can  be  written  as  equations  (9) 
and  (90  according  to  the  chemical  reaction  rate  laws: 


n 


T  d(flpt— p) 


Pt  d  t 


=  &4fCp0pt  — 


^4b^Pt-P  +  &5fCp0Pt-H  exp  ^  5 


-^5bCH+^pt— p  exp  - 


(1  -a5)n5Fr]^ 
RT 


+  kafCp0p{-u2  -  ^6b^Pt-pCn2  -  &7f#Pt-p 
(  uiqFrja 


x  exp 


x  exp 


(°n qFrh\  ,  y  n 

(  RT  J  +  k7bCH+0Pt_P^ 

(1  -  a7)qFq a\  (  a%qFq^\ 

- — I exp  I-—) 


+  k^CpC^+Opt  exp  (— 


(1  -  aQqFrh 
RT 


(9) 


d(0pt— p' )  7  a  ( aT‘l^rl/\  1  r'cl  a 

r Pt  df  =  ^7f#Pt-P  exP  I  RT  J  -  &7bCH+0pt-p' 


x  exp 


(1  -  cn)qFri 


a 


RT 


(90 


Substituting  equations  (6)-(8)  into  equation  (9)  for  0P{,  0pt_H 
and  0p{-n2  will  result  in  an  equation  (10)  which  only  contains 


0pt_p  and  ^pt-p'  as  the  surface  coverage: 


F 


T  d(flpt— p) 


Pt  dt 


=  A 


&4fCp  +  exp 


(1  -  a8 )qFq 


a 


RT 


+Bk6(CP  +  Ck5(CP  exp  h 


-<  A 


LtfCp  +  AxbC’p-C^*  exp 


-(1  -  ag )qFrj 


a 


RT 


+Bk6iCp  +  Ck5fCP  exp  +  k4 h 

(  —(l—aOnsFria  \ 

+^5bCH+  exp(  - — -  +^6bCH2  +%  exp 


RT 


( a%qFrh\  ,  .  / <xiqFr\a 

x  (  RT  J  +  h(  exp  I  —  )  }  0pi-p 


-  I  A 


&4fCp  +  A'xbC’p'  C‘4  exp 


-(1  -as,)qFq 


a 


RT 


+Bk6(Cp  +  Ck5(Cp  exp  Pppf  ~  kibC‘f  exp 


-(1  -  ai)qFr\a 
x  1  - ^ -  )  (  ^Pt-P' 


(10) 


The  differential  equations  (10)  and  (9')  can  be  solved  with 
the  boundary  conditions  of  that  at  t  =  0,  0P{-P  and  0pt_p'  =0, 
and  Cp  and  Cp'  =  0.  The  contaminant  surface  coverage  can  be 
obtained  as  a  function  of  contaminant  concentration,  anode 
overpotential  and  the  lifetime,  as  expressed  by  equations  (11) 
and  (IF): 


0pt-p 


F\ 


F\  F 2  +  F\  B7  +  B\ B2 


x 


1  (  F1B2+B1B2  —  B2B2+2F1F2—F2B2 

b^~2  I  ~~ 


F\  +  B\  +  F2  +  B2  —  Z 
x  exp  1 - - - 1 


1  /  F\  B2  +  B\  B2  —  B2B2  +  2F\  F2  —  F2B2  \ 

—Ab2+  -  j 


Fi  +  Bi  +  F2  +  B2  +  Z 
X  exp  | - r - 1 


(11) 


0pt-p' 


Fi 


F\  F2  +  F\B2  +  B\ B2 


x 


„  1  1  f  n  F\  B2+B1B2  —  B2B2  +  2F1F2  —  F2B2 

f2+~  (  b2- 
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X 


x  exp 


F\  +  B\  +  F2  —  B2  +  Z 
W\  -  B2) 

F\  +  B\  +  F2  +  B2 


-  Z 


1  /  F1B2  +  B1B2  —  B2B2  +  2F\F2  —  F2B2 

H —  Bi  H - 

2  V  Z 


x 


F\  +  B\  +  F2  —  B2  —  Z 


x  exp 


2(F\  -  B2) 

Fi  +  Bi  +  F2  +  B2  + Z 


t 


an 


where 


F\=  A 


ktfC?  +  kzbCP'Cl+  exp 


-(1  -  aq)qFi} 
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RT 


+BkbfCp 

B\  =  fc4b  +  k5bCH+  exp  —  j  +  ^6bCH2 

(  a7qFqa\  (  asqFqa 

+k-jf  exp  — — —  +  %  exp 


Fi  =  kj(  exp 


V  RT 
(  ajqFrk 
V  RT 


V  RT 


-(1  -  a7)qFr) 
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RT 


B2  =  k7bC{[+  exp 
Z  =  y/(Fi  +  ¥2  +  By  -  B2)2  -  4F2(Fi  -  B2) 


(11") 


In  equation  (11"),  F\  represents  the  sum  of  the  rate  con¬ 
stants  of  the  forward  direction  for  the  reactions  (IV)-(VI). 
These  forward  direction  reactions  are  mainly  responsible  for 
the  production  of  Pt-P.  B 1  represents  the  sum  of  the  reac¬ 
tion  constants  whose  corresponding  reactions  ((IV)-(VI)  and 
(VIII))  will  result  in  the  decrease  in  the  surface  coverage  of 
Pt-P.  F2  is  the  forward  rate  constant  of  the  reaction  (VII) 
which  is  responsible  for  the  increase  in  the  surface  coverage 
of  species  Pt-P7  and  partially  for  the  decrease  in  the  surface 
coverage  of  species  Pt-P.  B2  is  the  backward  rate  constant  of 
the  reaction  (VII)  which  is  the  only  contributor  to  the  decrease 
in  the  surface  coverage  of  species  Pt-P7.  Zin  (1 1)  and  (1 17)  is 
a  special  parameter  for  solving  the  differential  equations.  It 
can  be  proven  mathematically  that  Z’s  value  is  always  larger 
than  zero  if  F\ ,  B\ ,  F2  and  B2  all  have  non-zero  values  which 
is  the  case  proposed  above.  However,  if  F\  =  B2  (could  be  a 
possible  case  for  the  reaction  mechanism),  equation  (11)  will 
become  meaningless  (a  number  divided  by  zero).  If  this  case 
happens,  another  two  equations  (12)  and  (127)  rather  than 
equations  (11)  and  (1 17)  have  to  be  used  to  obtain  the  surface 
coverage: 


Opt-p  =  _  ,  p1  ,  „  [1  -  Fi  exp(-(Fi  +  Bi  +  F2)t)\ 

F\  +  B\  +  F2 

(12) 


Time,  Hour 

Fig.  2.  Calculated  surface  coverage  distribution  of  contaminant  P  and  its 
oxidation  product  P'  on  the  platinum  surface.  The  parameter  values  used 
are:  F\  =0.5,  B\  =0.2,  F2  =  0.1,  and  B2  =0.01  (the  unit  for  F2  or  B2 

is  mol  cm-2  h-1). 


0Pt-P'  — 


F2 


F\  +  B\  +  F2 
Fi 

x 


1  T-  _ 

B\  +  F2 

F\  +  B\  +  F2 
B\  +  F2 


exp(-(Fi  +  B\  +  F2)t) 


exp(-Fn) 


(12’) 


In  Appendix  A,  a  general  mathematical  discussion  about 
the  conditions  under  which  Z-  0  and/or  F\=B2  can  be 
found.  Although  several  cases  are  impossible  from  a  reac¬ 
tion  physics  standpoint,  the  discussion  is  still  necessary  for 
understanding. 

As  an  example,  Fig.  2  shows  the  calculated  distribution  of 
the  surface  coverage,  0pt-p,  ^pt-pr  and  0pt_p  +  0pt_p'  based  on 
parameter  values  ofFi=0.5,^i=0.2,  ^2  =  0.1,  and  B2  =  0.0 1 
(the  unit  for  F 1,  B 1,  F2  or  B2  is  molcm_2h_1).  It  can  be 
seen  that  at  the  early  stage,  both  0pt_p  and  ^pt-F  increases 
with  time.  After  a  while,  0pt-p  will  gradually  drop  and  0pt_F 
will  continue  to  increase  till  both  reach  steady- state  levels. 
The  drop  of  0pt_p  and  the  continuous  rising  of  ^pt-p'  reflect 
the  surface  electrochemical  conversion  from  species  Pt-P  to 
Pt-P7.  The  magnitude  of  the  steady-state  levels  is  determined 
by  the  magnitude  of  the  back  reaction  constant  (#1  and  B2)  of 
the  reactions.  When  the  time  is  long  enough,  the  surface  will 
be  largely  covered  by  Pt-P  and  Pt-P7,  that  is,  the  sum  of  0pt-p 
and  0pt-F  is  close  to  1 .  The  remaining  surface  available  for  H2 
electrochemical  oxidation  would  be  only  1  —  0pt-p  —  0pt-F , 
which  is  a  very  small  portion  of  the  whole  anode  surface.  The 
magnitude  of  1  —  0pt-p  —  0pt-F  is  also  mainly  determined  by 
the  value  of  the  sum  of  B\  +B2. 

In  the  case  that  P7  produced  by  the  electrochemical  oxi¬ 
dation  of  contaminant  P  is  not  adsorbed  (or  has  very  week 
adsorption)  on  the  platinum  surface,  the  surface  coverage  of 
species  P7  will  be  equal  to  zero  and  the  corresponding  items 
related  to  reaction  (VII)  will  disappear  from  the  equation  (9). 
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Solving  the  only  differential  equation  (9)  will  result  in  a  sim¬ 
pler  expression  for  surface  coverage  0pt-p: 

0Pt_P  =  [1  -  exp(-(F,  +  BOOl  (13) 

Fi  +  Bi 

A  typical  example  for  this  case  may  be  the  contaminant 
CO  in  the  hydrogen  stream  where  the  adsorption  of  product 
CO2  on  the  platinum  is  negligible  [25-30]. 

For  anode  H2  oxidation  (H2  — >  2H+  +  2e~),  as  discussed 
above,  reaction  (II)  may  be  considered  as  the  rate-determining 
step.  The  current  density  (4)  for  reaction  (II)  on  the  avail¬ 
able  catalyst  Pt  surface  could  be  expressed  as  equation  (14) 
[25,26]: 


4  =  ^H^Xa(fef^Pt-H2^Pt  -  feb^Pt-H) 


(14) 


Combining  equations  (6)-(8)  and  (14),  the  expression  for  fuel 
cell  anode  current  density  can  be  obtained: 


4  =  nuFy^f 


hf 

hb 


(14’) 


Table  1 


Parameters  and  their  corresponding  values  for  fuel  cell  performance  calcu¬ 
lation  in  the  absence  and  presence  of  anode  stream  contaminant  at  80  °C  and 
3  atm 


Parameter 

Value 

an 

2.0 

F 

96487  C  mol"1 

Ka 

30000  cm2  cm-2 

R 

8.314  JKmol-1 

T 

353  K 

no 

4.0 

Ydl 

5.1  x  10-4  Acm-2 

6^0^q'o 

0.5 

v° 

1.17V 

R° 

0.1  ^  cm2 

hfCn2/k\b 

0.1 

n0  /n0 
t7Pt/t7Pt— H 

2.0 

hf 

4.4  x  10-9  mol  cm-2  h-1 

hf 

5.0  x  10-2  molcm~2  ppm-1  h_1 

^4b 

6.4  x  10-2  mol  cm-2  h~ 1 

hf 

4.0  x  10-1  mol  cm-2  ppm-1  h-1 

^5b 

1.0  x  10-2  cmh-1 

hf 

1.0  x  10-5  molcm~2  ppm-1  h_1 

^6b  Cu2 

1.0  x  10-5  mol  cm-2  h-1 

hf 

0.0  mol  cm-2  h-1  ’ 

hb 

0.0mol1_<7  cm3<?-2h-1 ' 

hf 

8.0  x  10-8  mol  cm-2  h_1 

hb 

0.0  mol_<?  cm3<?+1  h- 1 1 

«3,5,7-8 

0.5 

n5 

1.0 

4 

2.0 

Assuming  the  corresponding  reaction’s  contribution  is  negligible. 


gen  reduction  (O2  +  4e  +  4H+  2H2O): 


where  4  is  the  fuel  cell  anode  current  density  and  ya  is  the 
ratio  between  the  real  Pt  surface  and  the  geometric  surface 
area  (cm2  Pt  cm-2  electrode)  for  the  fuel  cell  anode  catalyst 
layer.  For  a  fuel  cell  electrode  reaction,  the  reaction  zone 
is  a  three-dimensional  porous  catalyst  layer.  If  we  assume 
that  the  anode  hydrogen  oxidation  is  totally  controlled  by 
the  kinetics  (interface  diffusion  of  reactants  not  considered), 
a  simple  approach  is  to  introduce  a  parameter  to  reflect  the 
contribution  of  the  real  catalyst  reaction  area  to  the  current 
density.  Ao  is  the  value  of  A  expressed  by  equation  (6)  at 
the  equilibrium  electrode  potential  (^a  =  0),  which  could  be 
expressed  as: 


It  is  worthwhile  to  mention  here  that  for  a  fuel  cell  polar¬ 
ization,  a  flooded-agglomerate  model  has  been  employed  to 
describe  the  relationship  between  the  current  density  and  the 
overpotential  [32-34].  The  detail  of  this  model  will  not  be 
pursued  in  this  paper. 

For  the  cathode  side  in  the  absence  of  the  contaminant,  the 
relationship  between  the  current  density  (4)  and  the  cathode 
potential  (r/c)  can  be  given  by  equation  (15)  for  fuel  cell  oxy- 


rdo 


exp 


f  uonao  Ft]c 
V  RT 


—  exp 


f  ~(1  -  oto)naoF'qcY 

V  ft  )_ 

(15) 


Given  that  the  overpotential  for  oxygen  reduction  is  rela¬ 
tively  large  (>60  mV),  a  Tafel  equation  can  be  employed  to 
describe  the  cathode  polarization: 


0c  = 


RT 

P'OLq  F 


[ln(/c)  -  ln(yc*o)] 


(16) 


where  the  electron  transfer  coefficient  ao  and  the  electron 
transfer  number  nao  should  be  those  for  the  rate-determining 
step  of  oxygen  reduction  mechanism,  /q  is  the  exchange  cur¬ 
rent  density  of  oxygen  reduction  and  the  yc  is  the  ratio  be¬ 
tween  the  real  Pt  surface  and  the  geometric  surface  area  for 
the  cathode  Pt  catalyst  layer.  The  cathode  overpotential,  ijc,  is 
defined  as  the  difference  between  the  cathode  potential  (Ec) 
at  4  7^  0  and  the  equilibrium  potential  ( E [?)  at  4  =  0,  that  is, 
0c  ~  E°c  —  Ec. 

The  relationship  between  the  current  density  and  the  over¬ 
potential  has  been  calculated  according  to  the  parameter 
values  listed  in  Table  1  for  the  case  in  the  absence  of  con¬ 
taminant  (0pt-p  +  0pt-p'  =  0).  Kinetic  constants  were  chosen 
which  gave  a  reasonable  fit  and  trend  to  experimental  data 
reported  in  the  literature.  More  detailed  experimental  results 
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Fig.  3.  Calculated  fuel  cell  performance  and  corresponding  voltage  losses  of 
the  anode,  cathode,  and  membrane  according  to  the  parameter  values  listed 
in  Table  1  using  equations  (14),  (16)  and  (17)  (80  °C,  3.0  atm). 

for  the  contamination  of  the  anode  catalyst  will  be  acquired 
in  the  next  phase  of  research.  In  order  to  obtain  the  overall 
fuel  cell  performance,  another  parameter,  Rq  the  internal  re¬ 
sistance  of  the  cell  is  also  listed  in  Table  1.  Ignoring  mass 
transport  effects,  the  overall  cell  performance  can  be  calcu¬ 
lated  according  to  equation  (17): 

Tcell  —  V  7? a  0c  ^O^cell  (17) 

where  Vceii  is  the  fuel  cell  voltage,  V°  the  open  circuit  volt¬ 
age  and  /Ceii  is  the  current  density.  The  calculated  overall 
fuel  cell  performance  is  shown  in  Fig.  3  together  with  the 
anode,  and  cathode  polarizations,  and  ohmic  loss.  This  com¬ 
pares  well  with  other  performance  related  data  in  the  literature 
[10,45]. 

Substituting  equations  (11)  and  (IF)  into  equation  (14) 
for  0pt-p  and  0pt_ps  the  anode  overpotential  as  a  function 
of  time,  current  density,  and  contaminant  concentration  in 
the  presence  of  contaminant  can  be  obtained.  If  combining 
equations  (14),  (16),  and  (17),  (11)  and  (IF),  the  transient 
fuel  cell  voltage  in  the  presence  of  contaminant  can  be  sim¬ 
ulated  by  adjusting  reaction  constants  %,  k /b,  0pt/^pt-H  an<^ 
probably  ai  if  the  current  density,  anode  overpotential  and 
contaminant  concentration  are  experimentally  known.  This 
simulation  based  on  the  measured  cell  voltage  transient  be¬ 
havior  can  also  provide  a  method  to  estimate  the  chemical 
and  electrochemical  kinetic  constants. 

Note  that  the  resultant  transient  anode  overpotential  (r/a) 
and  cell  voltage  (VCell)  expressions  are  the  implicit  functions 
of  anode  overpotential.  It  is  difficult  to  get  an  analytical  ex¬ 
pression  for  ria  of  Vcell  •  Therefore,  some  numerical  calcula¬ 
tions  are  necessary  in  order  to  obtain  the  transient  behavior. 
As  examples,  for  a  simple  case  where  0pt_p'  =  0,  Fig.  4  shows 
the  calculated  results  at  various  stream  contaminant  concen¬ 
trations  based  on  equation  (13)  and  Fig.  5  for  those  at  various 
current  densities.  The  parameters  used  for  Figs.  4  and  5  are 
those  listed  in  Table  1 . 


Fig.  4.  Calculated  fuel  cell  voltage  degradation  in  the  presence  of  various 
levels  of  the  fuel  stream  contaminant  at  80  °C,  3.0  atm  and  a  current  density 
of  0.5  A  cm-2.  The  reaction  rate  constants  and  other  related  parameters  used 
for  this  calculation  are  those  listed  in  Table  1 . 


From  the  calculated  Figs.  4  and  5,  it  can  be  seen  that 
the  cell  performance  will  decrease  in  the  presence  of  an¬ 
ode  feed  stream  contaminant  until  a  steady- state  perfor¬ 
mance  is  reached.  The  greater  the  concentration  of  the  fuel 
stream  contaminant,  the  greater  the  rate  of  decrease,  and 
loss  in  performance  will  be,  and  the  sooner  steady-state 
performance  is  reached.  The  effect  of  current  density  is 
similar  to  contaminant  concentration.  The  fuel  cell  voltage 
degradation  trends  in  the  presence  of  fuel  stream  contami¬ 
nants  compare  well  with  other  reported  data  in  the  literature 
[27,29,30,46].  Again,  more  detailed  experimental  results  for 
voltage  degradation  will  be  acquired  in  the  next  phase  of 
research. 


Time,  Hour 

Fig.  5.  Calculated  fuel  cell  voltage  degradation  at  various  current  densities 
in  the  presence  of  5.0 ppm  contaminant  in  the  fuel  stream  at  80  °C  and 
3  atm.  The  reaction  rate  constants  and  other  related  parameters  used  for  this 
calculation  are  those  listed  in  Table  1 . 
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Current  Density,  A.cm’2 

Fig.  6.  Calculated  fuel  cell  steady- state  performance  as  a  function  of  fuel 
contaminant  concentration  and  current  density  at  80  °C  and  3.0  atm.  The 
parameter  values  used  for  the  calculation  are  those  listed  in  Table  1 . 


2.4.  Steady -state  performance  in  the  presence  of 
contaminant 

Figs.  4  and  5  shows  that  there  exists  a  steady-state  perfor¬ 
mance  for  each  contaminant  concentration  and  current  den¬ 
sity  when  the  lifetime  goes  to  infinity.  A  steady  anode  overpo¬ 
tential  and  performance  can  be  solved  numerically  according 
to  the  expression  of  steady-state  surface  coverage  obtained 
through  equations  (11)  and  (llOatt^oo: 


(1  -  0pt-P  -  0Pt— P')f=O0  = 


Ff  Ff  +  Ff  Bf  +  B*BS 


08) 


1  2 


1  2 


where  Ff,  Ff,  Bf  and  Bf  are  the  previous  defined  parame- 
ters  but  at  steady-state  under  the  condition  of  t  — >  oo.  Fig.  6 
shows  the  calculated  steady-state  polarization  of  a  fuel  cell 
at  various  contaminant  levels.  The  parameters  used  for  Fig.  6 
are  those  listed  in  Table  1 .  It  can  be  observed  from  Fig.  6  that 
the  steady-state  performance  level  is  a  function  of  contami¬ 
nant  concentration  and  current  density  (or  overpotential)  as 
mentioned  above. 

Note  that  if  the  rate  constant  of  reaction  (VII)  is  very  larger 
(larger  Ff)  or  the  sum  of  the  backward  rate  constant  is  negligi¬ 
ble  (B\  =  0),  there  will  be  no  steady- state  performance  which 
can  be  observed  due  to  the  entire  catalyst  surface  being  to¬ 
tally  occupied  by  P  +  P;  (that  is,  1  —  0pt-p  —  0pt_p'  ~  0).  In  this 
case,  the  cell  voltage  and  the  current  density  will  both  drop 
to  zero. 

For  a  special  case  where  the  reaction  (VII)  does  not  exist, 
the  equation  (18)  will  become  equation  (19): 


(1  -  0pt_P  -  0Pt_pO/=oo 


B 


s 

1 


Ff  +  Bf 


A  performance  loss  parameter  (PL%)  in  the  presence  of 
contaminant  can  be  defined  as  equation  (20): 


(20) 


Where  Vcsell  is  the  steady- state  cell  perfomance,  and  is 

the  steady-state  cell  performance  in  the  absence  of  contami¬ 
nation. 


2.5.  Contamination  transient  time  constant 

A  transient  time  constant  for  the  contamination  process 
can  be  defined  as  r,  which  is  a  measure  of  how  long  it  takes 
the  performance  to  approximately  reach  its  steady- state  level. 
The  expression  for  x  could  be  defined  by  equation  (21): 

r  = - - -  (21) 

F\  +  B\  +  F2  +  #2 

It  can  be  observed  from  equations  (11),  (IF)  that  the  tran¬ 
sient  time  constant  is  a  function  of  contaminant  level  and  an¬ 
ode  overpotential  (or  current  density).  In  general,  the  higher 
the  contaminant  concentration  or  current  density,  the  shorter 
the  transient  time. 


2.6.  Performance  recovery  process 


It  has  been  observed  experimentally  that  during  the  con¬ 
tamination  process,  the  fuel  cell  performance  could  be  recov¬ 
ered  automatically  with  time  when  the  contaminant  source 
was  cut  off  [2,3].  One  assumes  that  the  contaminant  source 
is  switched  off  at  a  time  of  to,  the  boundary  condition  for 
solving  equation  (10)  and  (90  would  be:  when  t  >  to,  Cp  =  0. 
The  obtained  solutions  are  given  in  equations  (22)  and  (220 
which  are  only  valid  under  the  condition  that  t>to'. 


i 


Ff- F2  +  F\-u<B2  +  BxB2 


X 


F[-t0B2- 


a 


x  exp 


2Zf-f° 

Ff'o  +  Bi  +  F2  +  B2-  Z^o 


(t  -  to) 


b 


2Zt-‘o 


+  Bi  +  F2  +  B2  +  Z‘^‘0 
x  exp  | - - - - - (t  -  to) 


(22) 


at>t0 

tfpt-p' 


F^>to 


2(F\-U]  -  B2)(F[-t0F2  +  F[^B2  +  B\  B2) 


(19) 


x 


F[~t(>  F2(F{-U]  -  B2) 
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+ 


F{- **>  +  Bi  +  F2  -  B2  +  Z1^ 


2Z,-,° 


F^o  +  Bx+F2  +  B2- 

x  exp  ( - L- - - - (t  -  to) 


+  (  F{^o  +  Bi  +  F2-B2-  Z^t(>  \ 


2Z‘^o 


F™  +  Bi  +  F2  +  B2  +  Z‘^ 
x  exp  I - - - - - (t-  to) 

(22’) 


where 

a  =  [(9p7-p (F[-t0  +  Bi  +  F2-B2-  Z'-t0) 

+2&K %'(F[-t0  ~  B2)\(F[-t0F2  +  F[-t0B2  +  BiB2) 
-2 F{~t0 F2(F[-t0  -  B2)  -  F{-t0  B2(F[-t0 
+F2  +  Bi  -  B2  -  Z'-r°) 


Fig.  7.  Calculated  surface  coverage  (0Pt_p  +  0Pt_P/ )  change  before  and  after 
to  at  various  F\ ,  B\ ,  F2,  B2  and  to  values  for  the  purpose  of  comparison.  The 
units  for  those  parameters  are  mol  cm-2  h_1 ,  except  for  to  which  is  in  hour. 
For  the  recovery  process  ( t  >  to),  F\  is  set  to  0  for  all  cases. 


b  = 


— 0  +  Bi  +  F2  -  B2  +  Z 
+2et=f,(F[-t(>  -  B2)\(F[-t0F2  +  F{-t0B2  +  BiB2) 
+2 F{-t0 F2(F{-‘°  -  B2) 

+F{-!oB2(F[-0  +  F2  +  Bi-  B2  +  zf-f0) 


T7t>tn  AVr  ~  nq  /-(l  -as)qFr]a\ 

F\-  =  AKsbCp'C^  exp  I  - — -  I 

Z‘-t0  =  y/iFf1 0  +  F2  +  Bi  -  Bi)-  -  4 F2(F{^  -  B2) 

where  and  0^^  are  those  surface  coverage  at  t- 
which  can  be  expressed  by  equations  (11)  and  (11')  except 
that  the  time  variable,  t ,  in  these  two  equations  is  replaced  by 
to.  F[-t0  is  the  F\  value  at  Cp  =  0.  If  the  product  concentration 
(CpO  produced  by  the  contaminant  oxidation  is  negligible, 
F\-t{)  will  be  approximately  equal  to  zero,  and  equations  (22) 
and  (22')  will  become  simpler.  Zt>‘o 

is  the  Z  value  at  Cp  =  0. 
The  similar  mathematical  discussions  to  the  equations  (11) 
and  (IT)  for  equations  (22)  and  (22')  are  listed  in  Appendix 
A. 

For  the  special  case  that  reaction  (VII)  does  not  exist, 
which  corresponds  to  the  case  where  the  surface  coverage 
of  P;  is  equal  to  zero,  the  surface  coverage  of  P  can  be  ex¬ 
pressed  by  equation  (23)  at  t  >  to’. 

0p?-p  =  f<><!+  Yi  {Fl?~  °  -  lptF°  -  CM- 0  +  *1)] 

x  expC-CFp'o  +  W-fo))}  (23) 

In  the  case  that  Cp  =  0,  which  corresponds  to  that  Fj-'0  =  0, 
equation  (23)  simplifies  to  equation  (24): 

$Pt— p  =  $pt— p  exp(— Bi(t  —  to))  (24) 


Fig.  7  shows  some  calculated  cases  with  various  forward 
and  backward  reaction  rate  constants  for  the  purpose  of  com¬ 
parison. 

The  results  show  that  the  larger  the  reaction  constant  for 
the  desorption  of  P  on  the  surface  (larger  B 1),  the  faster  the 
recovery  rate  would  be  (curve  1  versus  curve  2).  A  larger 
reaction  rate  constant  for  the  production  of  surface  P;  (larger 
F2)  would  make  the  recovery  slower  (curve  1  versus  curve 
3).  For  a  larger  B2  (the  rate  constant  for  the  surface  electro¬ 
reduction  of  P'),  the  recovery  process  would  become  faster 
(curve  3  versus  curve  4).  Another  observation  from  Fig.  7 
is  that  the  longer  the  contaminant  exposure  (longer  to),  the 
recovery  process  (curve  3  versus  curve  5). 

Combining  equations  (22)  and  (22')  with  equations 
(14)— (16)  allows  the  fuel  cell  voltage  to  be  solved  numerically 
for  the  recovery  process.  Fig.  8  shows  the  fuel  cell  voltage 
recovery  at  different  contamination  levels,  to’s,  and  current 
densities  for  the  purpose  of  comparison.  In  the  calculation, 
the  contribution  from  reaction  (VII)  has  been  omitted.  It  is 


Time,  Hour 

Fig.  8.  Calculated  fuel  cell  voltage  recovery  process  at  different  levels  of 
contaminant,  current  densities,  and  exposure  times  at  80  °C  and  3.0  atm. 
The  parameter  values  used  for  the  calculation  are  those  listed  in  Table  1. 
When  t  >  to,  F 1  is  set  to  zero  for  all  cases. 
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obvious  from  Fig.  8  that  Cp,  to,  and  current  density  can  make  a 
difference  for  the  fuel  cell  voltage  recovery  process.  A  higher 
contaminant  concentration  (1  versus  2),  a  higher  current  den¬ 
sity  (2  versus  3),  or  a  longer  exposure  (3  versus  4)  make  the 
recovery  slower. 

2.7.  Temperature  dependence  of  the  contamination 
process 

In  general,  temperature  can  affect  the  contamination  pro¬ 
cess  through  its  impact  on  the  parameter  values  in  the  equa¬ 
tions  discussed  previously.  An  increase  in  the  temperature 
will  cause  a  change  in  the  magnitude  of  the  reaction  con¬ 
stants  (kif  and  k ;b)  and  other  parameters  such  as  solubility 
of  hydrogen  and  contaminant  in  the  electrolyte,  their  diffu¬ 
sion  coefficients,  hydrogen  redox  exchange  current  density, 
electrolyte  resistance,  and  hydrogen  surface  recovery  in  the 
equilibrium  state.  For  the  contamination  adsorption  reactions 
(IV)-(VII),  rising  temperature  will  increase  the  magnitudes 
of  the  reaction  rate  constants  in  both  directions.  However,  the 
backward  reaction  rate  constants  would  be  increased  more 
than  that  of  forward  reaction  rate  constants,  resulting  in  a 
slower  contamination  rate  at  higher  temperature  than  that  at 
lower  temperature.  At  this  stage,  the  quantitative  approach 
for  the  temperature  effect  on  the  contamination  process  will 
not  be  pursued  in  this  paper. 


be  the  main  cause  of  the  anode  overpotential  increase  through 
the  elimination  of  the  active  Pt  catalyst  sites. 

The  fuel  cell  current  density  expression  as  a  function  of 
anode  and  cathode  overpotential  is  derived  from  the  proposed 
reaction  mechanism.  The  fuel  cell  performance  degradation 
in  the  presence  of  feed  stream  contaminant(s)  is  formulated 
as  a  function  of  contaminant  concentration,  current  density, 
and  lifetime.  The  obtained  equations  can  be  used  to  simulate 
and  estimate  the  chemical  and  electrochemical  reaction  rate 
constants,  and  make  some  prediction  about  the  severity  of  the 
contamination  and  the  performance  recoverability. 

Further  work  will  be  focused  on  modeling  of  the  mem¬ 
brane  and  cathode  poisoning  processes  and  validating  the 
models  with  more  detailed  experimental  data. 
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Appendix  A 

A.l.  t<to 


3.  Conclusion 

The  model  proposed  in  this  paper  for  the  anode  poison¬ 
ing  process  stems  from  the  reaction  mechanism  for  hydro¬ 
gen  oxidation  on  the  platinum  surface.  In  the  mechanism  of 
hydrogen  oxidation,  the  dissociative  chemical  adsorption  of 
adsorbed  H2  is  considered  as  the  rate-determining  step  for 
H2  oxidation.  The  electrochemical  reaction  of  dissociated 
H2  following  this  rate-determining  step  is  believed  to  have 
a  Nernst  behaviour  from  which  the  surface  coverage  of  H2, 
atomic  H  and  unoccupied  Pt  sites  are  derived  as  a  function  of 
contaminant  surface  coverage.  Chemical  adsorption  of  con¬ 
taminants  on  the  platinum  surface  and  their  impact  on  the 
electrochemical  hydrogen  surface  reaction  are  considered  to 


The  differential  equations  can  be  written  as  in  equation 
(A.l)  for  the  surface  coverage: 

d^pt-p 


d  t 

dOyt-p' 
d  t 


—  ~(F  1  +  B\  +  C2)0pt-p  —  (F\  —  Z^^pt-p'  +  F\ 


=  (F^pt-p  -  Z?2)0pt— p' 


(A.l) 


Regarding  the  initial  conditions  of  that  0pt-p  and  0pt-p'  both 
equals  to  zero  at  t  =  0,  the  analytical  solutions  for  equation 
(A.l)  can  be  obtained.  Five  cases  have  to  be  considered  in 
order  to  avoid  the  solutions  mathematically  meaningless. 

Case  1.  F\—B2^  0,  and  ( F\  +  B\  +  F2  —  Bf)1  +  47*2 

(F 1  —  B2)  7^0.  The  analytical  solutions  are  as  in  equation 
(A.2): 


0pt-p= 


F\ 


F\  F2  +  F\  Z?2  +  B\ Z?2 


B2- 


1  (  F\B2  +  B\B2  —  B2B2  +  2F\F2  —  F2B2  \ 

B2 - - -  )  exp 


F\  +  B\  +  F2  +  B2  —  Z 
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F\B2~\~  B\B2  —  B2B2  +  2F1F2  —  F2B2  \ 

B2  + - - -  )  exp 


Fi  +  Bi  +  F2  +  B2  +  Z 


0pt-p'  — 
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F\  F2  +  F\B2  +  B\B2 
-  exp 


1  /  F\  B2  +  B\ B2  —  B2B2  +  2F\  F2  —  F2B2  \  f  F\  +  B\  +  F2  —  B2  +  Z 

F2  +  -  (  b2 - - -  )  ' 


2(Fi  -  B2) 


F\  +  B\  +  F2  +  B2  —  Z  \  1 

- 5 - >  +rl«2  + 


F1B2  +  B\B2  —  B2B2  +  2FiF2  —  F2B2  \ 
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2 (Fl  ~  B2) 


exp 


F1  +  B1  +  F2  +  B2  +  Z 
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where  Z={{F\+Bi+F2-B2)2 +AF2{F\-B2))m  A-2-  *>*(> 


Case  2.  F\  —  B2  i=-  0,  B\  =  0  and  F\  =  F2  +  B2.  The  solutions 
are  as  in  equation  (A. 3): 

1 

0pt-p  =  —[ft  +  F2  -  (F\  +  Fi  -  F\  Fit)  exp(-fV)] 

F\ 

6»Pt_p'  =  ^[1  -  (1  +  Fxt)  exp(-fV)] 

F\ 

(A.3) 


Case  3.  F\  —  7^  0,  F2  =  0  and  B2  =  F\  +  B\ .  The  solutions 

are  as  in  equation  (A. 4): 


The  differential  equations  are  the  same  as  those  expressed 
by  equation  (A.l).  However,  the  boundary  conditions  for 
solving  the  differential  equations  are  different.  The  bound¬ 
ary  condition  would  be:  if  t  >  to,  Cp  =  0  and  0pt-p  =  0pf_5p 
and  0Pt_P/  =  0*^0, . 

The  analytical  solutions  for  equation  (A.  1)  can  be  obtained 
and  five  cases  have  to  be  considered  in  order  to  avoid  the 
solution  mathematical  meaningless.  Note  that  the  obtained 
equations  listed  below  are  only  valid  when  the  time  is  longer 
than  to- 

Case  1.  F\—B2^  0,  and  (F\  +  B\  +  F2  —  B2)2  +  4F2 

(F 1  —  B2)  7^  0.  The  analytical  solutions  are  as  in  equation 
(A. 7): 
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^p^_p/  —  0 


Case  4.  F\  —  B2  =  0  and  B\  +  F2  7^  0.  The  solutions  are  as 
in  equation  (A.5): 

Fi 
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F2  +  B\ 


exp(— Fit)] 


Case  5.  Fi  —  B2  =  0  and  B\  +  F2  =  0.  The  solutions  are  as  in 
equation  (A. 6): 
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(A. 6) 
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0p7_p  and  0^0,  are  those  surface  coverage  at  t-  to,  which  can 
be  expressed  by  equations  (18)  and  (180  except  that  the  time 
f  s  in  these  two  equations  are  replaced  by  Jo’s.  F[-t0  is  the  F\ 
value  at  Cp  =  0  and  Zr-r°  is  the  Z  value  at  Cp  =  0. 


Case  2.  F\  —  B2  7^  0,  B\  =  0  and  F\  =  F2  +  B2.  The  solutions 
are  as  in  equation  (A. 8): 
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(A.8) 


Case  3.  F\  —  B2  i=-  0,  F2  =  0  and  B2  =  F\  +  B\.  The  solutions 
are  as  in  equation  (A. 9): 

^-p  =  1  -  pt>to  ,  R  +  exP  (-Ff°(t  ~  to)) 

P\  + 

(  _  F^o  \ 

+  K4-  -  ffv  -  -pd  +  B)  exp  (A.9) 

x  (-(F[-f0  +  BOO  ~  to)) 
tfpjlp/  =  Of/ p,  exp (-F[-t0(t  -  to)) 


Case  4.  F\  —  B2  =  0  and  B\  +  F2  ^  0.  The  solutions  are  as 
in  equation  (A.  10): 
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+  1  +  B.  +  F, 


exp(— F[-t0(t  -  t0))  - 


F2 


B\  +  F2 


(A.  10) 


x 


rf — ^0  _ 

Ppt-p 


F 


t>to 

1 


exp(-(F[-r°  +  Bi  +  F2)(t  -  f0)) 


f‘-°  +  b{  +  f2 


1 


Case  5.  F\  —  B2  =  0  and  B\  +  F2  =  0.  The  solutions  are  as  in 
equation  (A.  11): 

#p7-p  =  1  +  (#p7-p  -  !)  exp(— -  t0)) 

(A.ll) 

0f%,  =  0 exp(-Fj-,0(f  -  f0)) 
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